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We propose a laser-driven Sisyphus cooling scheme for atoms confined in a far off resonance optical 
dipole trap. Utilizing the differential trap-induced ac Stark shift, two electronic levels of the atom 
are resonantly coupled by a cooling laser preferentially near the trap bottom. After absorption of 
a cooling photon, the atom loses energy by climbing the steeper potential, and then spontaneously 
decays preferentially away from the trap bottom. The proposed method is particularly suited to 
cooling alkaline-earth-like atoms where two-level systems with narrow electronic transitions are 
present. Numerical simulations for the cases of 88 Sr and 174 Yb demonstrate the expected recoil and 
Doppler temperature limits. The method requires a relatively small number of scattered photons 
and can potentially lead to phase space densities approaching quantum degeneracy in sub-second 
timescales. 

PACS numbers: 



I. INTRODUCTION 

Ultracold atomic gases are now routinely used for a 
broad range of science including precision measurements 
[lHil, quantum degeneracy studies @-[6j|, and quantum 
information applications 0, The standard proce- 
dure to produce such gases relies critically on laser 
cooling where the fundamental energy loss unit corre- 
sponds to the momentum recoil from photon absorption. 
The magneto-optical trap (MOT) arrangement for laser 
cooling can produce temperatures of hundreds of mi- 
croKelvins or even lower. Here, the combination of weak 
confinement and density dependent losses from photon- 
assisted collisions and radiation trapping restricts atomic 
densities to less than about 10 11 cm -3 [9(. For applica- 
tions requiring higher phase-space densities, evaporative 
cooling in magnetic or optical dipole traps (ODT) is used. 
While successful, evaporative cooling suffers from inher- 
ent atom loss and demands long trap lifetimes usually on 
the order of a minute. Most measurements with ultracold 
atoms can benefit from larger atom numbers and phase- 
space densities, as well as reduced experimental cycling 
times. 

In this paper we propose a method of Sisyphus cool- 
ing of atoms confined in an ODT that exploits dif- 
ferential ac Stark shifts of two electronic states. Our 
method promises phase-space densities which are much 
higher than those typically achievable with conventional 
laser cooling and in timescales which are much shorter 
than those typical of evaporative cooling. In particular, 
our scheme is well suited for rapid cooling of optically 
trapped spin-zero atoms which possess narrow electronic 
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intercombination transitions. 

In Sisyphus cooling the fundamental energy loss unit is 
the energy difference between two coupled atomic states 
at the positions of excitation and de-excitation, and can 
in principle be very large. Under suitable conditions, the 
atom is forced to repeatedly "climb uphill" and thus lose 
kinetic energy. The first such proposal by D. Pritchard 
[To| involved two rf-coupled hyperfine states in a mag- 
netic trap. Since then this idea has been investigated for 
a number of atom trap configurations. 

The original Pritchard scheme was theoretically stud- 
ied for the case of Ioffe-Pritchard traps [ll[ and modified 
to incorporate different magnetic energies of two elec- 
tronic levels fl2fl. Gravitational Sisyphus cooling in a 
magnetic trap [13[ and Sis yph us cooling in a blue de- 
tuned evanescent wave trap [14j have been demonstrated. 
RF-induced Sisyphus cooling in an ODT was demon- 
strated [l5j and a variant which exploits the second or- 
der Zecman effect and spin-exchange collisions has been 
proposed Recently a three-level cooling scheme ap- 
plicable to magnetically trapped Hydrogen has also been 
proposed [I?} • All of these schemes rely on transitions 
between internal magnetic sub-levels of the atom. None 
of these schemes is applicable to atoms with spin-zero 
ground state such as strontium or ytterbium. 

The rest of this paper is organized as follows. In sec- 
tion II we discuss the general requirements for our cool- 
ing scheme and determine suitable ranges of ODT wave- 
lengths for the application of our scheme to Sr and Yb 
atoms. In section III we develop a simple theoretical 
model to determine the efficiency of our cooling method 
and discuss its limitations. In section IV we develop a 
numerical model that properly addresses the stochastic 
nature of the process. We use this model to compute the 
expected value of final temperatures for Sr and Yb. In 
section V we discuss possible applications of our cooling 
scheme and draw our concluding remarks in section VI. 
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FIG. 1: (color online) Cooling scheme for a two-state atom 
in an ODT. An additional cooling beam, near resonant with 
the atomic transition, is also directed at the atoms. The solid 
upper red (lower blue) curve shows the spatially-varying ac 
Stark shift potential U e (U a ) of the excited (ground) atomic 
state induced by the trapping beam. The excited state polar- 
izability is greater than that of the ground state. The dashed 
blue curve is the ground state potential shifted up by the cool- 
ing photon energy. The frequency of the cooling beam is set 
to be resonant with the atom at ODT bottom. The cooling 
proceeds via the following cycle: 1. a moving atom absorbs 
a photon at the bottom of the ODT; 2. it climbs a steeper 
potential (solid upper red curve); 3. it spontaneously decays 
into the ground state with the shallower potential (solid lower 
blue curve). The atom loses energy e due to this process; 4. 
after a number of oscillations, the atom absorbs another pho- 
ton near the bottom of the ODT and the cycle repeats. 



II. COOLING SCHEME 

Our basic cooling scheme is shown in Fig. [T] A two- 
state atom confined in the optical dipole potential of the 
trapping laser interacts with cooling laser light that res- 
onantly couples the internal states near the trap bottom. 
A moving atom initially in the ground-state absorbs a 
cooling photon at the bottom of the ODT. The atom in 
the excited state climbs a steeper potential, thus spend- 
ing additional energy. Then it spontaneously decays into 
the ground state far from the minimum of the potential. 
After a number of oscillations the atom absorbs another 
cooling photon and the cycle repeats. 

Far off resonant optical dipole traps are essentially con- 
servative and atoms will scatter photons only from the 
cooling laser beam. Elastic collisions between atoms in 
the ground state are not required for our cooling method 
but are also not harmful unless the atomic sample is in 
the hydrodynamic regime. Even though the dynamics 
of our proposed cooling process are very different from 
that of a conventional MOT, the two share the same 
density limiting mechanisms - inelastic light-assisted col- 
lisions [lH and reabsorption of scattered photons. Cru- 
cially, we expect larger equilibrium densities than that 
achievable in MOTs due to the small excited state frac- 



tion as well as the fast cooling timescale. Together with 
the stronger spatial confinement, the scheme potentially 
allows far greater phase space densities than those achiev- 
able by conventional laser cooling. 

A. General requirements 

In order to implement the proposed scheme, one needs 
to satisfy the following conditions: (1) the ODT wave- 
length has to provide a higher polarizability for the atom 
in the excited state than in the ground state, (2) ab- 
sorption of the cooling light should occur only near the 
trap bottom, and (3) atoms in the excited state should 
have enough time to move substantially far from the min- 
imum of the ODT potential (ideally until turning point) 
to make the cooling process efficient, i.e. the natural 
decay time should be comparable to the trap oscillation 
period for the excited state. 

We consider the common case of loading of an optical 
trap directly from a MOT. Condition (1) can be satisfied 
by proper choice of the ODT wavelength. Condition (2) 
places a stronger restriction on the ODT wavelength, re- 
quiring that the relative ac Stark shift of the two states 
arising from their unequal polarizabilities is much larger 
than the linewidth of the cooling transition. This allows 
the cooling laser to be tuned to the atomic resonance only 
at particular locations in the trap, such as the bottom. 
Condition (3) implies an ideal natural decay rate com- 
parable to typical ODT trapping frequencies of less than 
100 kHz. For a particular choice of atom, conditions (1) 
and (2) may be satisfied with an appropriate choice of 
ODT wavelength. Condition (3) however mandates the 
availability of a narrow cooling transition. 

Our scheme is therefore well-matched to the case of 
alkaline-earth-like atoms (such as Ca, Sr, Yb) that pos- 
sess closed, spin-forbidden intercombination transitions 
with narrow linewidths. Ultracold samples of such atoms 
are of great interest in the context of optical clocks JlSj- 
l2l| . precise tests of fundamental physics 0, HH, |23j . 
quantum computing [24], [25[ and quantum simulation 
[26}. We focus on applying our scheme to the cases of 
88 Sr and 174 Yb and perform relevant calculations for the 
1 Sq — > 3 Pi transition (wavelength As r (Yb) = 689(556)nm, 
linewidth 7s r (Yb) = 2tt x 7.4(182) kHz). Laser cooling of 
these atoms on this transition is well developed. We first 
identify suitable ODT wavelength ranges that meet the 
requirements of our proposed scheme. 



B. Suitable ODT for Sr and Yb 

To determine ODT wavelengths that satisfy the above 
criteria for Sr and Yb, we calculate the ac Stark shift 
(see Fig. [2]) of the 1 5 , and 3 Pi levels for a large range of 
frequencies detuned to the red of the strong 1 So —> 1 Pi 
transition. The calculation of the polarizability is per- 
formed by summing up the light shift contributions with 
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FIG. 2: (color online) ac Stark shift of the 1 So (blue solid line) 
and 3 Pi (red dashed line) levels versus wavelength for 88 Sr and 
174 Yb. For concreteness, the calculation is performed for the 
maximum intensity in a single 1 W beam focused to a waist 
of 10 /xm. 



electronic states up to n = 8 principal quantum number 
for both Sr and Yb, using the spectroscopic data pre- 
sented in [53] . We have neglected the effect of hyperpo- 
larizability and considered a linearly polarized ODT. For 
the excited state we assume that atoms are in magnetic 
sublevels m = ±1 (28|. The range 700-800 nm appears to 
be suitable for Sr, since the polarizability of the excited 
state is substantially higher than that of the ground state. 
The ranges 690-800 nm and 1700-2500 nm are suitable 
for Yb. In these regimes, the differential ac Stark shift 
can be made much larger than the natural linewidth for 
either atom using modest ODT powers. Hence the cool- 
ing beam can be frequency tuned to address atoms only 
at the minimum of the ODT potential. 



III. ANALYTICAL MODEL 

To gain intuition about the dynamics of the pro- 
posed cooling process and understand its limitations, we 
first develop a ID analytic model where we neglect the 
Doppler shift and assume photon absorption strictly at 
the bottom of the ODT (x = 0). These assumptions will 
be relaxed in subsequent discussions. 



In our simple picture, we first estimate the cooling ef- 
fect from a single photon scattering event. If an atom 
initially in the ground state absorbs a cooling photon 
at time t — and goes into the excited state, the 
probability to still be in the excited state at time t is 
Pit) = exp (—7*). Here 7 is the natural decay rate of the 
excited state or equivalently the linewidth of the cool- 
ing transition. The average energy reduction e due to a 
single scattering event within time td is then expressible 
as 



(U e (t) - U g (t))jeM-jt)dt, 



(1) 



where U e i g \ (t) is the potential energy in the excited 
(ground) state. 

Assuming a harmonic potential f7 e (g) = muJ 1(g) x2 
the atom's position before spontaneous emission varies 
as x(t) — [(2fc j BT)/(mWg)] 1 / 2 sin(u; e i). Here w e (p) is the 
trap frequency in the excited (ground) state, m is the 
atomic mass, and E = ksT is the initial energy (ki- 
netic + potential) of the atom. We will be interested in 
timescales that are much larger than 1/7. For td — > 00 
in Eq. [1] we obtain 



1 



2 1 + ( 7 /2w e ) : 



(2) 



where cx e r g \ is the polarizability of the excited (ground) 
state. This is the average energy reduction from scatter- 
ing a single cooling photon. Here a g /a e = ^g/cj 2 < 1. 

We can now draw simple conclusions from this inter- 
mediate result. The cooling efficiency, i.e. the energy re- 
duction per scattered photon, strongly benefits from high 
trap frequencies with e = 0.5ksT(l — a g /a e ) in the limit 
Lu e 3> 7. However, already at ui e — 7, we have an energy 
reduction of e = 0.4fcsT(l — a g /a e ). The energy reduc- 
tion per transition is proportional to fc^T, i.e. the full 
energy of the atom, in contrast to more conventional laser 
cooling schemes where the energy change per photon ab- 
sorption is determined by the photon momentum. The 
proposed cooling scheme demands only a small number 
of scattered photons, and thus allows the use of transi- 
tions with a moderate branching ratio. As will be shown 
later, substantial cooling can be achieved with few tens 
of photon scattering events within a fraction of a second. 

We now consider two intrinsic heating sources for the 
atoms. The first is caused by photon recoil which lim- 
its the final energy to ksT ICC = h 2 'k 2 /m 1 where k is 
the wavevector of the cooling laser light. The second 
is caused by the uncertainty in the position at which the 
atom gets excited. The probability of absorption versus 
position in the ODT can be written as 



p(r) = A- 



l + s + 



g-g s (r)-k-v 

7 /2 



(3) 



where A is the normalization constant, s — I / I sa ± is the 
saturation parameter (intensity in units of the satura- 
tion intensity) of the cooling laser, and 5 is its detuning 
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from the transition frequency of the free atom. <5s( r ) is 
the ODT induced differential ac Stark shift between the 
excited and ground states and v is the atomic velocity. 
Even ignoring the Doppler shift, it is impossible to drive 
atoms into the excited state strictly at the minimum of 
the ODT potential due to the non-zero linewidth of the 
transition. This leads to a temperature limitation similar 
to the Doppler limit for cooling in MOTs, unlike the po- 
larization gradient cooling technique that overcomes the 
Doppler limit. Neglecting the Doppler shift term and as- 
suming a homogeneous distribution of atoms in the trap 
we can perform a simple estimate of this heating effect. 
In ID, taking 6 — Ss(0) (i.e. the cooling laser is resonant 
at the bottom of the trap), the average energy gained 
during an absorption event is 



AE 



p(x) 



W > 2 



(4) 



s which evaluates to AE = yl + s(hr//2) = /csTbop, 
exactly corresponding to the Doppler temperature Toop- 
This can serve as an estimate of achievable temperatures. 
We can then expect the cooling to cease at either T rec or 
Tbop, whichever is higher. 

We now include the Doppler shift due to atomic motion 
to evaluate the spatial dependence of transition probabil- 
ities due to the cooling laser. FigJ3] shows representative 
examples of photon absorption and subsequent emission 
probabilities as functions of position for various energies 
for Sr. The main consequence of the Doppler effect is to 
shift the position of peak absorption away from the trap 
center (see Fig. [3Ja)). This is because while being sub- 
stantially smaller than the ac Stark shift of the atomic 
levels, the Doppler shift can easily be larger than the 
linewidth of the cooling transition ( 1 S'o — > 3 -Pi). Even 
though the position shift is small compared to the trap 
size, the cooling efficiency will be reduced. We also note 
that the effect is far less pronounced in Yb because of its 
wider cooling transition. 

A possible strategy to overcome the effect of the 
Doppler shift is to introduce an additional detuning of 
the cooling light (see Fig.^b)). In this way the cool- 
ing light can be resonant for the atoms near the bottom 
of the trap. Further refinements include implementing a 
combination of two or more cooling beams at different 
frequencies or a frequency chirp of the cooling beam to 
follow the changing Doppler shift as the atom is being 
cooled. 

In the limit of trap frequencies uj e much higher than 
spontaneous decay rates 7 (w e 3> 7) the emission proba- 
bilities will be peaked near the turning points of motion, 
thus much further from the trap bottom than the absorp- 
tion probabilities and cooling can be very efficient. When 
uj e -C 7 the profile of the emission probability is essen- 
tially the same as that of the absorption probability and 
efficient cooling is impossible. When ui e is comparable to 
7 (as in Fig. [3]), the emission probability is still greatest 
near the motional turning points and efficient cooling is 
possible. This is the regime of focus in this paper. 
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FIG. 3: (color online) Absorption (solid lines) and subse- 
quent emission (dashed lines) probabilities of cooling photons 
per unit length versus position in ODT (centered at position 
zero). The calculations are performed for Sr atoms moving 
along the radial dimension of an ODT formed at the 10 /J,m fo- 
cus of a 10 W beam at 750 nm. Curve normalizations take into 
account that the time spent by an atom at a certain position 
is inversely proportional to its velocity. The integral under 
each absorption curve directly corresponds to the probability 
to absorb a photon during one oscillation. Emission curves are 
proportional to the rate of emission. The curves correspond 
to atomic energies of 400 (red (outer) curves), 50 (green (mid- 
dle) curves), and 5 (blue (inner) curves) /jK and cooling laser 
detunings of (a) 5 = S s (0) and (b) S = <5 S (0) - 12 7 . 



While the above discussion has been restricted to a 
single trapped atom, we now consider inelastic collisions 
in a many-atom system. We can make a simple esti- 
mate of the density limitations of this effect by balancing 
the growth in density from the cooling and the loss in 
density from two-body inelastic collisions. Atom den- 
sity in a harmonic trap scales as T~ 3 / 2 . Thus the rate 
at which density grows due to the Sisyphus cooling is 
dn/dt = (dn/dT){dT/dt) = (-2>n/2T){l/k B )(dE/dt). 
For an energy loss rate proportional to the energy, we 
can take dE/dt = —E/t and r = (ksT) / '(fje) where 
/ is the average fraction of time spent in the excited 
state and e(oc ksT) is the average energy loss per photon 
scattering event. Hyperfme-changing collisions are non- 
existent for spin-zero atoms making light- assisted colli- 
sions the dominant inelastic process. This process leads 
to a density decay rate — 2/3/n 2 where (3 corresponds to 
the two-body loss rate when / = 1/2. Our j3 is analo- 
gous to that used to characterize light-assisted collisions 
in bright MOTs where / ~ 1/2 pi |29j. The density 
reaches equilibrium when the two rates are equal, giving 
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"oq = 3/(4t//3). In the simple situation described by 
EqU 7i oq = (3 7 /8/3)(l - a g /a e )/(l + h/2u e ) 2 ). We will 
make quantitative estimates of r for our scheme and the 
consequent n cq in Section V. 



IV. NUMERICAL MODEL 

In order to address the stochastic nature of the process 
and to include the effect of the Doppler shift we develop 
a semiclassical Monte-Carlo simulation. Our numerical 
model is applied to the dynamics of a single trapped atom 
and is based on the fact that for a small enough time 
increment, changes in the relevant parameters such as 
atomic position and velocity are negligible. We break 
atomic motion in the ODT into a set of discrete time- 
steps At, and at each step calculate the probability of the 
atom to make a transition (either absorption or sponta- 
neous emission). We calculate the new position and ve- 
locity assuming that the acceleration is constant during 
At. If a transition occurs during a particular step, the 
atom starts the next step in the other state and expe- 
riences the different trapping potential. At is chosen to 
be much smaller than all other timescales of the problem 
(trapping frequencies and spontaneous decay rate) and 
is 1/zs or less for all the calculations in this paper. We 
include the photon recoil but neglect atom-atom interac- 
tions in the simulation. 

For an atom in the ground state, we model the ab- 
sorption of a cooling photon as a random process that 
happens with probability given by the right-hand-side of 
Eq. |3] with A = (j/2)At. We set 6 = S s (0) + S D , where 
5d = brings the laser into resonance with a stationary 
atom at the trap bottom and a finite 5d can be used to 
mitigate the effect of the Doppler shift. We model the 
spontaneous decay of an atom in the excited state as a 
random process with probability defined by the sponta- 
neous decay rate. 



A. ID case 

We first employ this numerical model for Yb in a ID 
harmonic potential (see Fig. 0]). We use trap frequencies 
corresponding to the center of the transverse profile at 
the focus of a 750 nm ODT with 2 W power and 3 fim 
waist. In Fig. U{a) we present the evolution of the en- 
ergy of the trapped atom. The discrete steps up or down 
in energy correspond to events of absorption or emission 
of a photon. The exact shape of the curve and the final 
energy varies between different numerical iterations due 
to the stochastic nature of the process. In Fig. |U[b), we 
show the average of 100 iterations with the same initial 
conditions. The energy drops to 5.5 /iK within the first 
3 ms, comparable to Tbop = 4.3 /xK, and then stays es- 
sentially the same. Here and in the rest of the paper we 
will express the energy of the atom in /iK. 
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FIG. 4: (color online) Energy of Yb atoms versus cooling 
time. Here At = 0.1/is. Harmonic potentials are assumed 
with trap frequencies w 9 ( e ) = 2n x 55(74) kHz (see text). The 
cooling beam parameters are s = 0.5 and Sd = — 7- (a) shows 
the atomic energy during three individual runs (iterations) of 
the numerical model. Each curve shows the time evolution 
for a single initial energy. The sharp steps correspond to 
transitions from ground to excited state and vice versa, (b) 
shows the results for the same conditions with each curve 
averaged over 100 iterations. Average final energy is ~ 5.5 /iK 
in this example. 

Similar simulations can be performed for Sr (see Fig. 
[5]). The smaller linewidth of the cooling transition leads 
to lower Doppler temperature but also a longer cooling 
timescale. The harmonic potentials correspond to the 
transverse profile at the focus of a 750 nm ODT with 5 
W power and 10 /im waist. To mitigate the more se- 
vere Doppler effect due to the narrower linewidth of the 
cooling transition (see FigJ3|), we introduce 3 frequency 
components in the cooling beam. Energies of about 1 /iK 
are reached within 15 ms. 

We studied the dependence of final energy on detun- 
ing (Fig. [5]) for the Yb case, with an initial energy of 
500 /iK. The initial Doppler shift is about 3 7. Clearly, 
very small detunings will not be useful since they cannot 
mitigate this shift. Detunings that are much larger than 
the Doppler shift would also clearly not be useful. This 
would suggest optimum (red) detunings on the scale of 
7, which is verified by our simulations. For the Sr case, 
a clear study is impeded by the requirement of multiple 
frequency components for good cooling performance. 

We also investigated the dependence of the final en- 
ergy on the intensity of the cooling beam. High intensity 
causes saturation broadening that makes the position of 
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FIG. 5: (color online) Energy of Sr atoms versus cooling time. 
Each curve shows the evolution for a single initial energy and 
is averaged over 50 iterations. Here At = 1/is. Harmonic 



potentials are assumed with trap frequencies to 
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10(14) kHz (see text). The intensity of the cooling beam is s — 
10 equally split between frequency components corresponding 
to Sd = — {5, 15, 45}7. Average final energy is ~ 0.86 fiK in 
this example. 
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FIG. 6: (color online) Final energy of a Yb atom versus cool- 
ing beam detuning expressed as — fc/7. Here At = 0.1 /zs, 
s = 0.5, and initial energy is 500 pK. The final energy is eval- 
uated as the energy after 20 ms, averaged over 100 iterations. 
Each point is the mean and uncertainty of four such values. 



photon absorption less denned. On the other hand low 
intensity of the cooling beam makes events of photon ab- 
sorption rare, slowing down the cooling process. There- 
fore, it is meaningless to ask about optimal intensity of 
the cooling beam without specifying the interaction time. 
To study this dependence, we set the interaction time to 
5 ms. We show in Fig. [7] the final energy versus intensity 
of the cooling beam and compare it to the Doppler limit 



of vTTs(?J7/2). 



B. 2D and 3D cases 

After gaining some basic understanding we can at- 
tempt to apply our numerical method to higher dimen- 
sions. The present cooling method strictly requires that 
the atom pass through the minimum of the trapping po- 
tential. This will always be the case in ID. 2 and 3 
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FIG. 7: (color online) Final energy of a Yb atom versus the 
intensity (expressed as saturation parameter) of the cooling 
beam. Here At = 0.1 /is, Sd = —7, and initial energy is 500 
/jK. The final energy is evaluated as the energy after 5 ms, 
averaged over 20 iterations. Each point is the mean and un- 
certainty of four such values. Also shown is the Doppler limit 
of conventional laser cooling for a given saturation parameter 
(blue dashed line). 



dimensional cases however allow trajectories that never 
cross the minimum of the trapping potential. This will 
make efficient cooling impossible. We verified this in sim- 
ulations with isotropic 2D or 3D confinement, performed 
in polar or spherical coordinates, where we observed re- 
duction in the radial velocity component but not in the 
angular ones. In such a situation, mixing between the 
different degrees of freedom is desirable. This can be 
achieved either by means of collisions with other trapped 
atoms, or by anharmonic mixing. We concentrate our 
study on the second mechanism, since it does not imply 
any requirements on collisional properties of atoms or the 
density of the atomic sample. 

Anharmonic mixing is a process that couples atomic 
motion between different dimensions of a trap. This 
coupling enables redistribution of energy over the differ- 
ent dimensions. ODT potentials are intrinsically anhar- 
monic. Since atomic energy is not negligible compared to 
the trap depth, anharmonicity plays an important role. 
However mixing between radial and angular components 
does not occur in spherically symmetrical potentials. To 
introduce mixing, an elliptical trap shape is required in 
the 2D case. In the 3D case, an ODT consisting of two 
intersecting beams of different transverse sizes can pro- 
vide the needed mixing. One can expect longer cooling 
times as well as higher final energies for 3D traps. This 
is because of the presence of trajectories that do not pass 
through the center of the ODT, causing the atoms to 
spend less time at the minimum of the trapping poten- 
tial. This reduces the excitation rate as well as the aver- 
age energy lost per scattering event, thus lengthening the 
cooling process. Additionally, the probability for atoms 
to be driven into the excited state away from the po- 
tential minimum is higher than in the case of ID traps. 
This increases the average energy gained per excitation 
and limits the final energies to values higher than in ID 
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(see Eq. 0}. 

We consider two 2W ODT beams with wavelength 
750 nm, focused to waists of 3 and 3.5 fim. intersecting 
perpendicularly at their foci. Fig. [5] shows an example of 
cooling of Yb in such a trap. While a single cooling beam 
is sufficient to lower atomic energy substantially, our sim- 
ulations indicate that 2 or 3 beams [28[ will provide faster 
cooling and about 30% lower final energies. Employing 
3 orthogonal cooling beams at total intensity s = 0.5 al- 
lows energy decrease from 3 x 500 /iK to ~ 3 x 10 /iK 
during 40 ms of cooling from scattering as few as 20-40 
photons. As expected, 3D cooling takes a considerably 
longer time than in ID. We write the energy as lksT per 
dimension. This definition of T then directly corresponds 
to the temperature for a sample of many atoms. While 
thermal equilibrium is not guaranteed by our scheme, 
equilibrium will be established within a few collisional 
timescales after the cooling process. 




10 20 30 40 

time (ms) 



FIG. 8: (color online) 1/3 of total energy of Yb atoms versus 
cooling time, evaluated in 3D simulations. Each curve shows 
the evolution for a single initial energy and is averaged over 
50 iterations. Here At = 0.1 [is. The atoms move in an 
anharmonic crossed-optical dipole trap potential, as described 
in the text. The cooling beams have detuning 8d = —7 and 
total intensity s — 0.5, split into three equal parts aligned 
orthogonal to each other with two along the ODT beam axes. 
The average final energy after 50 ms of cooling is ~ 3 x 8.6 /iK. 

We summarize our findings in Table I where we report 
the final energies for Sr and Yb for ID and 3D simula- 
tions. The Sr energies (10xl2/im crossed ODT with 5W 
in each beam) were obtained for the saturation parame- 
ter s = 10 per cooling beam. Lower final energies could 
be obtained at the cost of longer cooling time at lower 
intensities. Our results compare well with the limits set 
by the Doppler and recoil temperatures. 





Jrec 


TDop 






Sr 


0.46 


0.18 


0.8 


3 x 1.4 


Yb 


0.35 


4.3 


6.0 


3 x 8.5 



TABLE I: Lhe final energies (in /iK) for Sr and Yb obtained 
in 1 and 3 dimensional simulations. Each value of final energy 
is the result of 80 time averaged numerical simulations. 



V. POSSIBLE APPLICATIONS 

Our proposed method can potentially enhance many 
experiments based on optically trapped atoms, includ- 
ing studies with quantum gases, precision metrology, 
and quantum information science. The method can act 
as an additional cooling step that improves initial, pre- 
evaporation conditions for these experiments. For prop- 
erly chosen parameters, it may be possible to increase the 
phase space density to approach quantum degeneracy. 
The achievable phase space densities can be estimated 
by incorporating the limitations from density-dependent 
loss processes. Optimum performance will require a care- 
ful choice of ODT and cooling beam parameters. It is 
experimentally easier to implement the proposed cool- 
ing method in ODTs with a larger volume. However a 
sufficiently high optical power at given wavelength is not 
always available. We perform calculations for the cases of 
174 Yb atoms in 2000 nm wavelength ODT and 88 Sr in 800 
nm wavelength ODT. High power at these wavelengths 
can be obtained from commercially available lasers. 

We first examine the case of 174 Yb atoms in an ODT 
formed by the perpendicular crossing of beams with 
waists 10 and 12 /im, each with 20 W power. For an 
initial energy of 3 x 500 /iK, our one-atom numerical sim- 
ulations indicate an average final energy of ~ 3 x 9 /iK 
using a single cooling beam with s — 3 and 8r> = — O.57. 
The trap depth (5mK) and volume seem adequate to ef- 
ficiently capture atoms from a MOT operating on the 
So — > 3 Pi transition. 

We can estimate the density limitations from light- 
assisted collisions using the expression n cq = 3/(4t//3) 
derived in Section III. Here we will assume that this is 
the dominant limiting effect. For simulations with initial 
energy 0.5 mK, we obtain average values of r = 65 ms 
and / = 0.0025. Using p = 10 _n cm- 3 / s [H, we then 
get n cq ~ 10 15 cm~ 3 as the limiting density. The re- 
sultant phase space density at a temperature of 9 /iK 
is 0.05, which greatly exceeds what is typically achiev- 
able in optically trapped Yb, prior to evaporative cool- 
ing [33|, [34|. While the estimated final optical density 
near the center of the trap is about 6, radiation trapping 
of the emitted photons should not be an issue, because 
emissions take place preferentially away from the center, 
at the motional turning points. Furthermore, the short 
timescale for cooling will alleviate losses from other in- 
elastic processes such as three-body collisions as well. For 
subsequent studies, the trapped sample lifetime can be 
extended by adiabatically relaxing the ODT confinement 
after cooling. 

The cooling performance can be improved by using 
time-dependent cooling beam parameters. Using 88 Sr 
atoms in an ODT formed by the perpendicular cross- 
ing of beams with waists 10 and 12 /jm, each with 10 
W power (7mK depth), our simulations indicate cooling 
from 3 x 250 /iK to ~ 3 x 1.15 /iK with three orthogonal 
cooling beams with parameters s — 0.5 + 5.0e _t / 100ms 
and 8d = — {3, 10, 20)7. These simulations indicate av- 
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erage values of r = 78.1 ms and / = 0.01, which yield 
n eq ~ 1.9 x 10 14 cm" 3 using (3 = lO" 11 cm" 3 /s HI 
and a consequent final phase space density of 0.5. This 
is of particular significance as direct evaporative cooling 
of 88 Sr is infeasible due to the near-zero atomic collision 
cross-section. The density limiting mechanisms may also 
be further mitigated by using time varying cooling beam 
detunings and relaxing the ODT confinement near the 
end of the cooling process. We also note that the sponta- 
neous scattering rate of ODT photons for the above cases 
are 0.3 s _1 (Yb) and 18s _1 (Sr). The large value for the 
Sr case may start to affect achievable temperature and 
phase-space density. 

We now turn to applications of this scheme beyond 
the wavelength ranges and atomic species discussed thus 
far. We first note that the use of time-dependent cooling 
beam parameters could allow extension of the proposed 
method to optical trap wavelengths satisfying a g /a e > 1. 
In this case, S has to be increased over time, starting from 
5-0. 

In addition to alkaline-earth-like atoms, this method 
can be employed for other atomic species as well, pro- 
vided sufficiently narrow cycling transitions are available. 
In alkali atoms, while the usual D2 line roSi/2 — > 
will be much too broad, the narrower nS\/2 — > (n + 
l)P 3 / 2 line may be utilized. This transition has recently 
been used for conventional laser cooling of 6 Li(3|| and 
4O K0. 

Application of this method to spin-polarized Fermi 
gases can prove a useful alternative to sympathetic cool- 
ing strategies. In this case an additional advantage is 
derived from the Pauli suppression of inelastic collisions 
between ground state atoms. A particularly appealing 
application may be in the context of an optical lattice 
trap where only a few atoms are located in each lattice 
site. Atoms in different sites will then enjoy high trap 
frequencies and thus efficient cooling, while at the same 
time the total number of traps and therefore atoms can 
be substantial. 



VI. CONCLUSIONS 

We have described a Sisyphus cooling method for 
atoms confined in an optical dipole trap that exploits 



the trap-induced differential ac Stark shift of two atomic 
levels coupled through a narrow-linewidth optical tran- 
sition. With a proper choice of ODT wavelength and 
energy levels, the atoms are cooled by a process of prefer- 
ential absorption of photons from a resonant laser beam 
near the trap center followed by spontaneous emission 
away from the center. The presented cooling scheme re- 
sembles other Sisyphus cooling methods but doesn't rely 
on the presence of ground state magnetic sub-structure 
and is well suited for cooling alkaline-earth-like atoms 
which have narrow intercombination transitions. 

Numerical simulations were presented for 88 Sr and 
174 Yb which show that temperatures as low as few /iK 
are reachable in timescales of tens of milliseconds. The 
temperature is limited to the higher of Doppler and re- 
coil temperatures. We pointed out experimentally acces- 
sible parameters for Yb and Sr where our predicted final 
phase space density from this cooling method is near the 
quantum degenerate regime. Further improvements in- 
clude the use of time dependent trap geometries as well 
as time dependent cooling beam parameters. The scheme 
is also adaptable to narrow electronic transitions in other 
atomic species. 

The presented cooling method falls under the more 
general class of one-photon cooling methods where no- 
ticeable reduction of the temperature is achieved by scat- 
tering of one photon, as demonstrated in [37], H|[ • These 
methods are crucial for cooling species that do not pos- 
sess near-cycling transitions with high branching ratio. 
An appealing prospect is cooling of molecules that have 
been pre-cooled using other methods [39h41| |. Recent 
work in the DeMille group has demonstrated laser cool- 
ing of heteronuclear molecules 42[ on a transition with 
a moderately high branching ratio, adequate for cooling 
using the proposed method. 
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